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We express the coverage (occupation fraction) 9, in powers 
of time t for four models of two-dimensional lattice random 
sequential adsorption (RSA) to very high orders by improving 
an algorithm developed by the present authors [J. Phys. A 
29, L177 (1996)]. Each of these series is, to the best of our 
knowledge, the longest at the present. We analyze the series 
and deduce accurate estimates for the jamming coverage of 
the models. 

PACS number(s): 05.50.+q, 05.20.-y, 05.70.Ln 



I. INTRODUCTION 

A large variety of physical, chemical, biological, and 
ecological processes can be modeled by random sequen- 
tial adsorption (RSA) |lj on a lattice. In this irreversible 
process, particles are deposited randomly on a surface 
one after another. In the simplest form, the depositing 
particles are represented by non-overlapping hard core 
extended objects. These particles, once adsorbed on a 
surface, are assumed to be fixed in their spatial positions 
and exclude certain regions from further occupation. It 
is therefore interesting to study the coverage (fraction of 
occupied sites) of the system as time evolves, in particu- 
lar the jamming coverage, a standstill beyond which no 
deposition is possible. 

Exact solutions for RSA models are usually restricted 
to special cases only H||. For those models which re- 
sist an exact theoretical treatment, different approximate 
methods have been developed |jj . Series expansion is one 
of the powerful analytical methods in RSA studies [p|-^0[ . 
In this work we present an improved algorithm for series 
expansions based on hierarchical rate equations P,|lO|] . It 
is then applied to four lattice RSA models where we are 
able to obtain, to the best of our knowledge, the longest 
series known for each of model that we have studied. 



II. MODELS 

Four models are studied in this work. For each of 
these models, we start with an initially empty, infinite 
lattice. Particles are dropped randomly and sequentially 



at a rate of k per lattice site per unit time, onto the 
lattice. Hereafter we set k equal to unity without loss 
of generality. In the first model, monomers are deposited 
on a square lattice with nearest-neighbor (NN) exclusion, 
where an adsorbed monomer excludes its occupied sites 
and its four nearest-neighbors from further occupation 
HPI. I n the second model, dimers of random orienta- 
tions are deposited on a square lattice. In this model, a 
dimer can occupy a pair of nearest-neighbor sites only if 
both of the chosen sites are vacant ||[llj . The third and 
the fourth models are defined similarly as the first and 
the second models, respectively, except that depositions 
are performed on a hexagonal lattice OJl3] • 



III. MODIFIED ALGORITHM 

We shall first summarize the original algorithm pre- 
sented in Ref ||. We are interested in expressing the 
coverage as a function of time t, which is related to 
the probability of finding a single vacant site (denoted 
by 'o'). More precisely, we write 9{t) = 1 — P(o,t). In 
general we are interested in expressing the probability 
of finding any particular configuration of interest, Co, 
which we write it as Pc = P(Co), as a power series ex- 
pansion of time t, with the expansion point t = 0. We 

write Pc (t) = S^Lo / n '' w ^ n derivative of 

d"Pc (t) 



P Co given by P^ 

Let Gi denote the set of new configurations generated 
in the calculation of the ith derivative of Pc , and G\ 
the corresponding jth derivatives of the set of configu- 



rations. We observe that Gq , G™" 
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mined at the (n - 1 )th derivative] , Gq~ 2 , G™~ 3 . 
[determined at the (n — 2)th derivative], . . ., Gq are pre- 
determined before calculating the nth derivative of Pc a ■ 
In the calculation of nth derivative of Pc a , we determine 
systematically Gq, G™~ , . . ., G\_ 1 , G° , by recursive use 
of rate equations. This algorithm is efficient since each 
value in G™ -1 , < i < n and the rate equation for each 
configuration G are generated once only. However, this 
algorithm consumes the memory quickly due to storage 
of intermediate results. In this algorithm, each configu- 
ration is transformed into a unique canonical represen- 
tation and stored as such in memory. If a configuration 
is already generated before, only a pointer reference to 
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the previous configuration is made. The checking of the 
existence of a particular configuration is done efficiently 
through the hashing technique ]l3| . 

An instructive way of looking at this algorithm is 
through the concept of tree traversal |lj| . We make the 
first configuration Co , of which its highest order of deriva- 
tive to be obtained is h, as the root node of the tree. The 
depth at which the root node resides is set arbitrarily to 
0, this corresponds to the index i = in the symbol G\ 
introduced above. We see that a complete breadth-first 
tree traversal will introduce at new set of leaves nodes of 
depth which is one more than the depth of leaves nodes 
in the previous breadth-first traversed tree. Therefore at 
the end of Dth times of tree traversal we will introduce a 
set of new leaves nodes of depth D. Since we have to store 
all new configurations or nodes that appear in the pro- 
cess, the breath-first tree traversal is thus a very memory 
demanding algorithm. To curb the exponential memory 
growth, we stop the breadth-first traversal of tree when 
the leaves nodes have a depth D' where a new and dif- 
ferent strategy is followed. Notice now that the problem 
of finding the hth order derivative of Co amounts to pre- 
determining the first, second, . . ., through (h — D') order 
derivatives for all the leaves nodes of depth D'. If we can 
do that, then by another (h — D') times of breadth-first 
traversal of the tree from the root node Co again, we can 
calculate all derivatives of Co up to order h. Now the 
problem is to calculate the first through (h — D')th order 
derivatives for each of the leaves nodes of depth D'. This 
can be achieved by treating each of the leaves nodes of 
depth D' as a new root node in another breath-first tree 
traversal subproblem. Each of these subproblems can be 
dealt with one at a time, which means that we can use 
the memory allocated for a subproblem for the next sub- 
problem. It is through this way that the use of memory 
is kept at roughly a fixed amount. A drawback of this 
approach is that the derivatives for some configurations 
have to be calculated more than once, however there is an 
important observation that the calculations of all these 
subproblems can be parallelized. 

As another effort to cut the growth of the memory re- 
quirement, we observe that the storing of the nodes is not 
really necessary for nodes which have an absolute depth 
very close to h. This can be achieved by using a depth- 
first tree traversal. When we are considering the rate 
equation of a particular configuration C, we recursively 
traverse down the tree whenever a configuration belongs 
to the right hand side of the rate equation for node C, 
which we shall call it a child node of C, appears. The 
merit of this depth-first approach is that we do not need 
to consider all symmetry operations that have to be ap- 
plied to a configuration. These operations are required 
before any configuration can be kept in the memory. Al- 
though the depth-first strategy effectively saves the trou- 
bles of finding a canonical representation, the time to 
traverse a depth-first tree can be quite long if the height 
of the tree is large. 

A detailed description of the implementation of the 



original algorithm can be found in Ref |10[. It can be 
modified readily to suit this new improved algorithm. 
Here we should mention the details of the task of per- 
forming the rate equation expansion. When the canon- 
ical representation for a new configuration is required, 
we make a list of 2-column vectors of entries which are 
the x and y coordinates of vacant sites specifying the 
configurations stored in a "working matrix" . Symmetry 
operations are applied to the new configuration through 
pre-multiplying the vectors in the list with appropriate 
transformation matrices. Criteria are set so that we can 
determine a canonical representation out of the list of 
new vectors. 

Series expansions are parallelized using PVM version 
3.3.11 Jl4| to utilize the hardware resources to the fullest. 
Calculations were performed on a 16- node cluster of Pen- 
tium Pro 200 with a bonded dual channel 100 BaseT 
Ethernet connection. All runs that we have done use less 
than 60 Megabytes of main memory, with the longest 
run took an elapsed time of 120 hours per machine in 
a dedicated computing environment. To compare the 
performance of this new algorithm and the original algo- 
rithm j^], we notice that we are able to obtain 4 more 
coefficients for the dimer RSA model on a square lattice 
using the former algorithm. The later algorithm requires 
706 Megabytes of memory even at the 14th order deriva- 
tive calculation. The results for all the models are dis- 
played in Table | and Table for square lattice RSA 
and hexagonal lattice RSA, respectively. We note that 
our 18th order coefficient in the expansion of powers of 
U = e~ t - 1, is 3560240252651011168 for monomer RSA 
with NN exclusion on a square lattice, which is slightly 
different from the value reported in Ref |sj . 

IV. ANALYSIS OF SERIES 

The approach to the jamming state for lattice RSA 
is often exponential. This behavior enables us to use 
transformations of variables which reflect the actual ap- 
proach. The transformations that we have employed are 
similar to that used in Refs (5^ . First we transform the 
coverage 9(t) in term of y = 1 — e . Another trans- 
formation z = (1 — e~ by )/b is performed. This transfor- 
mation is suggested by the fact that the exact jamming 
approach of dimer RSA on a linear lattice corresponds 
to 8{t) = z with b = 2 For other RSA models, 
we set b to be a free variable so that when different or- 
ders of Pade approximants are applied to 9(z), a cross- 
ing region between different orders of Pade approximants 
is to be located so as to give a good estimate of jam- 
ming coverage 0oo. As an example, Fig. ^ shows the 
crossing region of Pade approximants of orders [11,10], 
[10, 11], [12, 9], [9, 12], [13, 8], [8, 13], [14, 7], [7, 14], 
[10, 10], [11, 9], [9, 11], [12, 8], [8, 12], [10, 9], [9, 10], 
[11, 8], and [8, 11] for the series of monomer RSA with 
NN exclusion on a square lattice, giving an estimate of 
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6>oo = 0.3641323(1), where the last digit denotes the un- 
certainty. This estimate is in good agreement with the 
estimate of 6*00 = 0.3641330(5) by Baram and Fixman 
||. Analyzing the series using the square cactus as a 
reference model [|, we obtain 9^ = 0.364132(1). All 
these estimates from the series analysis agree well with 
the simulation result of 9oc = 0.36413(1) p5[ . 

For the dimer RSA on a square lattice, we obtain 9 X = 
0.906823(2). This is to be compared with the simulation 
results of Oliveira et al. § of X = 0.90677(6) and of 
Wang and Pandey |l6| of 9^ = 0.906820(2). A somewhat 
biased estimate of 9^ = 0.9068088(4) was obtained in 
Ref |)| as too small a range of b was used. 

For monomer RSA on a hexagonal lattice, 9^ = 
0.37913944(1). This result is much more accurate than 
the simulation results of 6™ = 0.38(1) and 9^ = 0.379 
reported in Refs JIt]] and Jl^] , respectively. For the dimer 
RSA on hexagonal lattice, we obtain doc = 0.8789329(1). 
This is in good agreement with the estimate 9 X — 
0.87889 in Ref Q. 

V. CONCLUSIONS 

In this work we have substantially extended the series 
for 9(t) by at least 3 terms for RSA models defined on a 
square lattice. The series for the hexagonal lattice RSA 
are also of very high orders. We noted that with the im- 
provement of the original algorithm, very long series and 
accurate estimates for 9^ can be obtained. The general- 
ity of this computational method allows us to handle a 
variety of problems based on rate equations. We are cur- 
rently applying this algorithm to provide more insights 
into the dynamics of the two-dimensional kinetic Ising 
model @. 
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TABLE I. Series-expansion coefficients for the probability of finding a vacant site of RSA models on the square lattice, with 
n designating the order of derivative. 
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l 


— 1 


-4 


2 


5 


28 


3 


-37 


-268 


4 


349 


3212 


5 


-3925 


-45868 


6 


50845 


756364 


7 


—742165 


— 14094572 


8 


12017245 


292140492 


9 


—213321717 


—6653993260 


10 


4113044061 


164952149516 


11 


-85493084853 


—4416119044972 


12 


1903886785277 


126863203272268 


13 


-45187885535477 


-3889473277203116 


14 


1137973688508989 


126677386324657804 


15 


-30289520203949205 


-4365431744153008620 


16 


849248887429012733 


158621692698054953164 


17 


- 250072598709248 1 7749 


-6058617368871081964076 


18 


771322713104711008093 


242593022530588935132428 


19 


-24860884250598911650645 




20 


835568036857675195155997 




21 


-29227671255970830546587445 





TABLE II. Series-expansion coefficients for the probability of finding a vacant site of RSA models on the hexagonal lattice, 
with n designating the order of derivative. 



n 


monomer RSA with NN exclusion 


dimer RSA 


1 


-1 


-3 


2 


4 


15 


3 


-22 


-99 


4 


154 


807 


5 


-1306 


-7803 


6 


12946 


87039 


7 


-146026 


-1097139 


8 


1837666 


15383607 


9 


-25429018 


-237018699 


10 


382667218 


3973893519 


11 


-6208467946 


-71934898755 


12 


107847476914 


1396766151303 


13 


-1994552336218 


-28932596959515 


14 


39089479606162 


636355373676831 


15 


-808620801445066 


-14801609728262739 


16 


17597107973050354 


362819105840203479 


17 


-401684435661166234 


-9343388660601611115 


18 


9593278364964882706 


252096661713936722415 


19 


-239168310954693706954 


-7109230409530804525155 


20 


6211745800773276251122 


209085655917886138504551 


21 


- 167765099068145498846842 


-6400562190944268227677947 


22 


4703750163644209363257538 


203577552403910264228987775 


23 


-136703709388456354551539146 




24 


4112464853135639607791658706 
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